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A REVIEW  OF  MAXIMUM-ENTROPY 
SPECTRAL  ANALYSIS 


R.W.  Herring 


ABSTRACT 

'J 

The  maximuTn-entropy  method  (MEM)  of  spectral 
analysis  is  a technique  for  estimating  the  power 
spectra  of  time-series  data.  The  purpose  of  this 
note  is  to  survey  the  literature  of  the  MEM  and 
to  provide  some  insight  into  its  derivation  and 
usage.  Fortran  subroutines  for  computing  the 
MEM  power  spectra  of  complex-valued  data  are 
included. 


1.  INTRODUCTION 

The  maximum-entropy  method  (MEM)  of  spectral  analysis  is  a technique 
for  estimating  the  power  spectra  of  time-series  data.  When  the  "true" 
spectrum  of  the  process  being  investigated  consists  of  discrete  lines 
separated  by  at  least  the  reciprocal  of  the  length  of  the  data  record,*  the 
MEM  yields  a spectral  estimate  consisting  of  sharp,  narrow  spikes.  This  is 
in  marked  contrast  to  the  performance  of  the  classical  Blackman  and  Tukey 
(1959)  and  periodogram  techniques.  These  latter  techniques  estimate  spectral 
lines  as  broad,  smooth  peaks  which  tend  to  merge  if  the  lines  are  closely 
spaced  in  frequency.  The  discontinuous  nature  of  the  spectrum  is  thus  masked. 


The  classical  Rayleigh  criterion. 
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The  purpose  of  this  report  is  to  survey  the  literature  of  the  MEM 
(Section  2)  and  to  provide  some  insight  into  its  derivation  and  usage 
(Sections  3 and  4).  Section  3 includes  a derivation  of  an  MEM  algorithm 
suitable  for  analyzing  complex-valued  data  records.  The  algorithm  is  not 
new  (cf.  Smylie  et  al.,  1973),  but  the  present  derivation  follows  the  more 
lucid  approach  taken  by  Andersen  (1974)  for  real-valued  data.  Section  3 also 
Includes  a brief  description  of  the  criterion  derived  by  Akaike  (1969a, b; 

1970)  for  selecting  the  most  appropriate  of  the  many  possible  choices  of  the 
MEM  spectrum  which  can  be  derived  from  a particular  set  of  data.  Section  4 
refers  to  listings  of  Fortran  subroutines  for  computing  the  MEM  spectra  of 
complex-valued  data.  These  listings  are  Included  as  appendices. 

MEM  spectral  analysis  of  complex-valued  data  is  applicable  to  data 
derived  from  such  systems  as  sampled-aperture  antenna  arrays,  where  the 
calculated  spectra  may  be  interpreted  in  terms  of  signal  strength  as  a function 
of  angle-of-arrlval.  The  method  should  also  be  useful  for  the  frequency 
analysis  of  tone-encoded  data.  In  such  applications  the  use  of  complex- 
valued data  allows  the  unambiguous  distinction  between  positive  and  negative 
frequencies,  which  often  is  necessary. 


If  -• 


2.  REVIEW 

The  two  classical  methods  of  spectral  analysis  are  the  periodogram 
technique  and  the  Blackman  and  Tukey  (1959)  power  spectrum  technique.  In 
using  the  periodogram  technique,  the  squared  magnitude  of  the  discrete 
Fourier  transform  (DFT)  of  the  data  is  computed.  Nowadays  the  DFT  is  usually 
computed  by  means  of  the  fast  Fourier  transform  (FFT)  algorithm.  The  ampli- 
tude (unsquared)  spectrum  is  often  smoothed  (usually  by  the  well-known 
technique  of  "windowing"  the  data  before  transformation)  in  order  to  reduce 
the  excessive  spreading  or  "leakage"  cf  spectral  lines  caused  by  the  finite 
length  of  the  data  record.  Often  the  computed  power  (squared-magnitude) 
spectrum  is  also  smoothed  to  Improve  its  statistical  reliability.  Both  types 
of  smoothing  reduce  spectral  resolution. 

The  Blackman  and  Tukey  technique  consists  of  two  stages.  In  the  first 
stage,  the  autocorrelation  function  of  the  data  is  estimated  by  extending  a 
finite-length  data  record  with  zeros  and  then  computing  the  autocorrelation 
function  of  this  extended  data  record.  Usually  lag  distances  of  not  more 
than  10%  of  the  original  record  length  are  considered.  In  the  second  stage 
of  the  procedure  the  DFT  of  the  estimated  autocorrelation  function  is 
computed  to  obtain  a smoothed  power-spectrum  estimate.  This  technique  was 
developed  before  the  FFT  algorithm  was  known,  when  the  direct  computation  of 
DFTs  was  very  laborious.  It  was  Intended  for  use  in  analyzing  noise-like 
data  with  smoothly  varying  spectra  (l.e.  , having  no  spectral  lines),  and  it 
avoided  the  costly  and  frustrating  necessity  of  first  computing  a spectrum 
in  far  greater  detail  than  was  required,  only  to  subsequently  smooth  it  in 
order  to  obtain  the  desired  statistical  stability. 

The  maximum-entropy  method  of  spectral  analysis  was  first  proposed  by 
Burg  (1967,  1968)  as  a technique  for  improving  spectral  resolution  in  the 
sense  of  reducing  the  spectral  smearing  caused  by  the  finite  length  of  the 
data  record.  In  the  classical  spectral  analysis  techniques  described  above. 
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the  estimated  autocorrelation  function  is  tacitly  extended  either  by  zeros  or 
by  a cyclic  repetition  of  Itself.  Burg  proposed  that  a more  appropriate 
extension  of  the  autocorrelation  function  could  be  based  on  that  extrapolation 
which  gives  rise  to  the  statistically  most  random  (i.e.,  maximum-entropy) 
power  spectrum  estimate  consistent  with  the  observed  data.  Burg  (1968)  showed 
how  to  calculate  such  a maximum-entropy  spectrum  from  a real  sampled  data 
record.  This  technique  and  its  extension  to  include  complex-valued  data  is 
presented  in  Section  3 which  follows. 

Van  de  Bos  (1971)  demonstrated  that  the  Burg  maximum-entropy  equations 
were  mathematically  the  same  as  those  encountered  in  the  least-squares  fitting 
of  an  all-pole  spectral  model  to  time-series  data.  He  also  pointed  out  that 
such  a model,  in  which  the  data  are  assumed  to  have  been  generated  by  a set 
of  damped  resonators,  is  invalid  for  many  kinds  of  data.  This  restriction 
must  be  taken  into  account  whenever  the  use  of  the  maximum-entropy  method 
is  proposed.  In  many  practical  cases*,  however,  the  data  can  be  considered 
as  the  sum  of  one  or  more  truncated  series  of  complex  exponentials  plus  noise, 
and  the  all-pole  model  is  appropriate. 

Lacoss  (1971)  did  a comprehensive  analysis  of  both  the  MEM  and  the 
analytically  somewhat  similar  maximum-likelihood  method  (MLM)  of  spectral 
analysis.  The  major  result  of  this  work  was  to  obtain  the  shape  and  magnitude 
of  the  expected  spectral-response  function  for  the  case  where  a small  number 
of  exactly  known  samples  of  the  autocorrelation  functions  of  real  sinusoidal 
signals  plus  uncorrelated  additive  noise  have  been  given.  The  crucial 
problem  of  estimating  the  autocorrelation  function  was  not  considered.  It 
was  also  shown  that  if  the  autocorrelation  function  were  known  exactly, 
sinusoids  with  frequencies  separated  by  amounts  greater  than  the  reciprocal 
of  the  length  of  sampled  autocorrelation  function  would  be  sharply  resolved. 
This  is  in  marked  contrast  to  the  blurred,  broad  spectral  peaks  which  occur 
in  this  case  for  the  classical  Blackman-Tukey  and  periodogram  methods.  It  is 
recommended  that  the  prospective  user  of  the  MEM  examine  the  results  of 
Lacoss'  work  closely,  particularly  in  regard  to  the  shape  and  magnitude  of 
the  expected  spectral  peaks. 

Ulrych  (1972)  and  Ulrych  et  al.  (1973)  devoted  considerable  effort 
to  the  problem  of  using  the  analysis  techniques  outlined  by  Burg  (1968)  to 
analyze  short,  real-data  records  containing  one  or  more  sinusoidal  components. 
Their  results  included  a comparison  of  the  Fourier  periodogram  and  maximum 
entropy  spectra  of  the  same  records,  and  showed  a dramatic  sharpening  of  the 
spectral  peaks  estimated  by  the  maximum  entropy  method  relative  to  those 
estimated  by  the  Fourier  method. 

Smylie  et  al.  (1973),  in  a very  complete  and  rigorous  exposition, 
generalized  Burg's  analysis  technique  to  include  the  case  of  complex  sampled 
data.  The  extended  technique  was  then  used  to  analyze  irregularities  in  the 
earth's  rotation.  An  incomplete  extension  to  complex  data,  based  on  the 
elegant  recursive  formulas  derived  by  Andersen  (1974) , was  published  by 
Haykln  and  Kesler  (1976) , who  seemed  unaware  of  the  earlier  and  more  compre- 
hensive work  of  Smylie  et  al.. 


e.g.,  short  records  of  tone-encoded  data;  plane-wave  data  received  by  a linear  sampled-aperture  array  of 
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The  most  comprehensive  review  to  date  of  the  subject  of  MEM  and  MLM 
spectral-analysis  techniques  Is  that  of  Ulrych  and  Bishop  (1975) , which  showed 
the  applicability  of  the  criterion  derived  by  Akalke  (1971)  for  determining 
from  the  data  the  optimal  length  of  the  spectral  estimator  (see  Section  3) . 

Another  excellent  review  of  the  subject  was  recently  published  by 
Kaveh  and  Cooper  (1976).  The  relation  between  the  MEM  and  the  very  similar 
autoregressive  method  Is  clearly  brought  out  (cf.  van  den  Bos,  1971)  and  a 
comparison  Is  made  between  these  methods  and  the  classical  Blackman  and  Tukey 
(1959)  approach. 

The  problem  of  estimating  confidence  Intervals  for  MEM  spectral 
estimates  has  been  considered  by  Baggeroer  (1976) , who  included  the  case  of 
estimating  the  confidence  limits  for  spectral  estimates  based  on  complex  data. 
This  highly  mathematical  paper  included  a comparison  of  the  results  of 
Akalke  (1969b) . 

Finally,  the  problem  of  estimating  maximum-entropy  spectra  when  a 
set  of  measured  autocorrelation  functions  of  unequal  statistical  weight  are 
available  has  been  Investigated  by  Newman  (1977).  This  extension  to  the  MEM 
is  particularly  useful  when  estimated  autocorrelation  functions  derived  from 
long  records  of  data  have  been  produced  by  special-purpose  digital  processors. 
In  this  case  a more  efficient  calculation  of  the  maximum-entropy  spectrum  can 
be  made  then  by  using  the  method  of  Burg. 


3.  THE  BURG  MAXIMUM-ENTROPY  METHOD* 

The  Burg  maximum-entropy  method  of  spectral  analysis  Is  particularly 
suited  to  the  analysis  of  short  data  records  when  estimates  of  the  auto- 
correlation function  have  not  been  given.  For  purposes  of  the  derivation, 
however,  it  is  easier  at  first  to  suppose  that  a series  of  M+1  equi-spaced 
complex  autocorrelation  function  values  4>(k),  k=0,l,...,M  are  known.  Here 
4>(k)  is  the  autocorrelation  with  lag  kAt  where  At  is  the  sampling  interval. 
From  the  usual  definition  of  the  autocorrelation  function,  4>(-k)  = (|)*(k), 
where  the  asterisk  denotes  complex  conjugation. 


The  power  spectrum  P(f)  estimated  by  the  maximum-entropy  method  is 
given  by 


P(f) 


"m"^ 

M 

1 - expf-j2x  fk  At) 

k=l 


2 


(1) 


i 


The  exposition  which  follows  is  an  extension  of  the  work  of  Andersen  (1974),  except  that  here  the  data 
are  assumed  to  be  complex  rather  than  real.  The  present  results  are  identical  to  those  obtained  by  the 
more  rigorous  derivation  by  Smylie  et  al.  (1973). 
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where  the  frequency  f is  limited  to  the  Nyquist  Interval,  -l/(2At)  £ f £ l/(2At), 
and  the  parameter  and  the  coefficients  are  determined  by  solving  the 
matrix  equation 


4>(0)  <j>*(l)  ...  (|)*(M) 

1 

r "m" 

<(>(1)  4)(0) 

"^1 

0 

: 4>*(i) 

• 

4>(M)  $(1)  <K0) 

0 

(2) 


Equations  (2)  show  that  Pj^  is  the  output  power  of  the  (M+l)-length  prediction- 
error  filter  (1,  -aj^i , . . . .-aj^) . An  (M+l)-length  prediction-error  filter 
gives  as  its  output  the  difference  between  the  datum  and  a linear  estimate 
of  its  value  based  on  the  proceeding  M data. 

The  concept  that  P(f)  as  given  by  (1)  is  a spectral  estimator  can  be 
understood  as  follows.  Equation  (2)  shows  that  the  output  of  the  (M+1)- 
length  prediction-error  filter  is  a single  pulse,  which  has  a white  power 
spectrum  Pn^t.  This  prediction-error  filter  is  thus  a filter  which  whitens 
the  power  spectrum  of  the  input  process.  The  reciprocal  of  the  power 
spectrum  of  the  whitening  filter  therefore  gives  an  estimate  of  the  power 
spectrum  of  the  process. 

In  general,  after  the  solution  of  (2)  for  order  M has  been  obtained, 
it  is  desired  to  extrapolate  to  order  M+1.  The  solution  for  order  M+1  in- 
volves the  determination  of  a set  of  M+3  unknowns  (aj^+i  i » • • • iSM+l ,M+1 > 

(})  (M+1)  ;Pm+i)  from  the  M+2  eqns . (2).  Thus  an  additional  criterion  is  required 
in  order  to  find  a unique  solution.  Burg  (1968)  suggested  that  this  criterion 
be  based  on  that  choice  of  which  minimizes  the  output  power  Pxi+1  > 

where  Pj^j+i  is  the  average  of  the  output  power  of  the  prediction  error  filter 
operating  on  the  data  in  both  the  forward  and  reverse  directions. 


For  the  trivial  case  M=0,  P^  is  estimated  by 


P 

o 


_1 

N 


N 

I 

t=l 


(3) 


where  P^  is  an  estimate  of  <Ji(0)  and  N is  the  number  of  complex  data  samples 
x^.  It  is  again  assumed  that  the  sampling  interval  is  At. 


For  the  case  M=l,  the  two-point  prediction-error  filter  (1,  -aj^i)  is 
found  by  minimizing  the  average  output  power  IT^  with  respect  to  a,,,  where 
is  given  by 


n 


1 


1 

2 N-1 


N-1 

I 

t = l 


11 


t+i' 


(4) 
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The  first  term  In  the  summation  Is  the  forward-direction  output  power  of  the 
filter  and  the  second  Is  the  reverse-direction  output  power.  The  use  of  the 
complex  conjugate  of  the  filter  coefficients  Is  a consequence  of  the  require- 
ment that  the  power  spectra  of  the  forward  and  reverse  direction  filters  be 
Identical.  If  the  data  were  real,  then  x*  = x , and  as  a consequence  the  a's 
would  be  real  as  well  (a*^^  = ^11^* 

Maximizing  as  a function  of  a^j^  yields 


a 


11 


N-1 

Z 

t=l 


N-1 

2 Z X*  X 


(5) 


This  step  completes  the  first  part  of  the  Inductive  proof  of  a recursive 
procedure  for  obtaining  the  sets  of  a’s  and  P's  from  a given  set  of  complex 
data.  Substitution  of  the  value  for  a^^^^  from  (5)  Into  eqn.  (4)  yields, 

* ^1’  '^h^ch  may  then  be  used  to  make  an  estimate  of  <(>(1)  from  eqns.  (2). 
In’orSer  to  evaluate  the  a's,  however,  explicit  estimates  of  the  (j>'s  are  not 
required.  This  fact  will  become  evident  below. 

A useful  set  of  recursion  relations  between  the  filter  coefficients 
aM-i  k k output  powers  Pm-1  ^nd  Pjj  will  now  be  derived,  follow- 

ing Burg  (196§).  From  eqn.  (2)  can  be  written  recursive  matrix  relations  to 
link  the  (M-1) ^^-order  filter  to  the  M*-^-order  filter,  when  M+1  values  of 
(J)(k)  are  known: 


<^(0) 


1 


^,M-1 

■®MM 


<l>(0)  

• • 

• • 

• • 

1 

"^-1,1 

-a 

MM 

0 

“®M-1,M-1 

•« 

• • 

• • 

• • 

• 

"®M-1,M-1 

-^-1,1 

((■(M)  

■'  'J'(O) 

0 

L J 

1 

(6) 


y 


: 


where 


\-l  = Vl,M-k 

k=l 

Equations  (6)  and  (7)  lead  immediately  to  the  recursion  equations 


^,k  ^=l,k  " ^^M-l,M-k’  ^ 1 ^ 1 


~ ^^M-1 


Additional  useful  definitions  are 


= -1 


= 0 (k  > M) 


From  eqn.  (9),  the  average  output  power  II^  of  the  M order  filter  is  then 


N-M  / M 


2 (N-M)  ^ I J'  ^Mk  ’^t+k  ®Mk  *t+M-k 

t=l  \ k=0  k=0 


N-M  / M 


2 (N-M)  I ^-l,k  ’'t+k  “ ^ ^-l,M-k  ^t+k 


M M 

®M-l,k  *t+M-k  " ^fM  , ^-l,M-k  ^'t4M-k 

k*0  ’ k*0  ’ 
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or 


1 / 
\ 2(N-M)  y 


where 


M M 

E a*  , . X.  .,  = E a* 


Hit  ° ^-l,k  *t+k  Hl-l,M-k  ^t4M-k 


'^Mt 


M 

E a. 


M 
= E 


1^=0  "t4W-k  Hl-l,M-k  ’'t+k 


t = 1,2 N-W 


Since  b„  and  b'  are  independent  of  a^,  the  condition 
Mt  Mt  MM 


Hi  = 


= 0 


gives 


N-M 


N-M 


Him 


= 2 E b*  b'  / E 

t=l  t=l 


(ivi'"  iSti') 


Furthermore 


a , N-M  / 

4 "«  ■ - efl  ('^ 


Mt' 


+ I’)  > 0 


so  that  the  extremum  of  for  a^^  given  by  eqns.  (14a)  and  (14b)  is  a 
minimum . 

Useful  recursion  formulas  for  the  arrays  b and  b'  are  derived 
means  of  (9)  and  (14a)  and  (14b) : 


'’Mt  “ 


I - a*  b' 

M-l,t  Hi-1,M-1  M-1,1 


b'  * b' 
Mt 


'M-l,t+l  ~ Hi-1,M-1  ^M-l.t+l 


(13) 

(14a) 

(14b) 

(15) 

(16) 
(17) 
by 

(18a) 

(18b) 


This  completes  the  second  part  of  the  inductive  proof.  It  has  now  been  shown 
that  the  filter  for  any  order  M can  be  constructed  from  the  filter  of  order 
M-1. 
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It  can  be  seen  that  the  arrays  and  bjl{t  are  constructed  from  the 
arrays  b^^_^  ^ and  ^ y t;  ^ simple  linear  operation.  The  starting  values 
are  ’ ’ 


"ot  = ’^Ot  = "t 

(19) 

but  in 
values 

practice 
are  used 

the  iterative  procedure  is  stated  at  M=1  and 
Instead : 

the  following 

b = X 

It  t 

(20a) 

b'  = X 

It  t+1 

(20b) 

t - 1,2, , ,N-1. 

From  eqn.  (16)  it  is  apparent  that  la^l  £ 1;  therefore,  from  eqn.  (9), 
it  follows  that  0 < P,  < , . 

The  Inverse  of  the  autocorrelation  matrix  4>  can  be  expressed  in  terms 
of  the  filter  coefficients  and  the  output  powers  (Burg,  1968).  If  ip  = 
then 


\li*  = = \b* 

^ji  ^M-i,M-j  ^ M-j,M-i 

n 

^-k,i-k  ' ^-k, j-k'^^M-k 

l,j  = 0,1 M, 


where  n is  the  smaller  of  the  two  Indices  i and  j . This  inverse  may  be 
required  in  certain  applications. 


A criterion  for  the  choice  of  the  appropriate  filter  length  M(<N)  has 
been  described  by  Akalke  (1969a,  b;  1970)  and  summarized  by  Ulrych  and  Bishop 
(1975) . The  criterion  is  known  as  teh  Akaike  final  prediction  error  (FPE) 
and  is  defined  as  the  mean-square  prediction  error.  Numerically,  the  FPE  of 
an  order  filter  is  given  by 


(FPE)^ 


N + (M+1) 

N - (M+1)  M 


if  the  sample  mean  has  been  eliminated  from  the  data,  or 


(22) 


(FPE)^ 


N + M 
N - M M 


(23) 
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otherwise*.  The  appropriate  choice  of  M Is  then  that  which  yields  the  minimum 
FPE.  There  appears  to  be  no  easy  way  of  determining  this  minimum  short  of 
direct  evaluation  for  each  value  of  M,  and  care  must  be  taken  to  avoid  local 
minima  Induced  by  statistical  fluctuations.  Most  authors  suggest  that  a 
restriction  of  M £ N/2  be  observed  to  prevent  the  use  of  excessively  long 
filters,  which  tend  to  artificially  split  what  ought  to  be  single  spectral 
peaks  (Chen  and  Stegen,  1974). t 


A final  observation  (Andersen,  1974)  Is  that  eqn.  (1)  gives  a detailed 
picture  of  the  whole  power  spectrum,  but  If  It  Is  desired  only  to  determine 
the  frequency  and  magnitude  of  the  principal  peaks  of  the  spectrum.  It  may  be 
more  appropriate  to  determine  the  positions  of  the  complex  roots  of  the 
pol)momlal 


P(z) 


(24) 


relative  to  the  unit  circle,  and  evaluate  (1)  only  In  those  regions  where 


4 


4.  COMPUTING  THE  MAXIMUM  ENTROPY  SPECTRUM 

Appendix  A contains  the  listing  of  a Fortran  subroutine  CMESA,  which 
computes  the  residual  powers  P2,P2»*-*»Pm  filter  coefficients 

k = 1,2,...,M  from  a given  set  of  complex  data. 

Appendix  B contains  the  listing  of  a Fortran  subroutine  CMESAP,  which 
computes  the  magnitude  of  the  power  spectrum  at  a particular  frequency 
according  to  (1) . 

The  roots  of  the  polynomial  p(z)  of  eqn.  (24)  can  be  calculated  using 
the  subroutine  CPOLY  (Jenkins  and  Traub,  1972;  Withers,  1974).  A tested  copy 
of  this  subroutine  package  Is  available  from  the  author  of  this  document. 


5.  CONCLUSION 

The  literature  of  the  maximum-entropy  method  of  spectral  analysis  has 
been  briefly  summarized,  and  the  more  significant  references  have  been 
singled  out  as  essential  reading  for  the  potential  user. 

The  Burg  method  (Section  3)  appears  to  offer  a way  of  resolving 
spectral  peaks  In  short  data  records  when  such  peaks  are  separated  by  at 
least  the  order  of  the  reciprocal  of  the  length  of  the  data  record.  Limited 
simulation  experimentation  by  the  author,  using  two  complex  sinusoidal  signals 


\ 


t 


* If  the  process  being  investigated  has  zero  mean  and  can  give  rise  to  signals  at  frequencies  near  zero,  then  the 
sample  mean  should  not  be  eliminated  from  the  data  and  eqn.  (23)  is  the  appropriate  formula  to  use. 

t But  see  also  Jones  ( 1 976). 


differing  In  frequency  by  about  a half  a cycle  over  the  record  length,  has 
indicated  that  for  such  data  the  computed  spectra  were  very  sensitive  to  the 
relative  phases  of  the  sinusoids.  In  particular,  good  resolution  was  obtained 
when  the  sinusoids  were  in  phase  quadrature  at  the  mid-point  of  the  record, 
and  spurious  results  were  obtained  in  all  other  cases.  This  result  apparently 
occurred  because  the  cross-product  terms  in  the  short  estimates  of  the  auto- 
correlation functions  did  not  cancel. 

The  development  of  the  algorithm  for  the  processing  of  the  complex- 
valued  data  has  been  given  in  detail  and  a Fortran  version  of  its  implementa- 
tion has  been  included.  It  is  hoped  that  this  Technical  Note  will  assist 
those  involved  in  the  evaluation  of  the  applicability  of  MEM  spectral  analysis 
to  their  particular  sets  of  data. 
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DO  310  T = 1,N-MM 

NUM  = NUM+C0NJG(B1(T) )»B2(T) 

DEN  = DEN+REAL(B1 (T) )#*2+AIMAG(Bl (T) )*#2 
+REAL(B2(T) )**2+AIMAG(B2<T) )*»2 
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